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Abstract 

The origin of the breaking of conventional hnear /c^-factorization for hard processes in a nuclear 
environment is by now well established. The realization of the nonlinear nuclear /c_l -factorization 
which emerges instead was found to change from one jet observable to another. Here we report on 
an important technical progress, the evaluation of the four-gluon color dipole cross section operator. 
It describes the coupled seven-channel non-Abelian intranuclear evolution of the four-gluon color- 
singlet states. An exact diagonalization of this seven-channel problem is possible for large number 
of colors Nc and allows a formulation of nonlinear A;_l -factorization for production of gluon-gluon 
dijets. The momentum spectra for dijets in all possible color representations are reported in the 
form of explicit quadratures in terms of the collective nuclear unintegrated glue. Our results fully 
corroborate the concept of universality classes. 
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I. INTRODUCTION 



The key point behind conventional perturbative Quantum Chromo Dynamics (pQCD) 
factorization theorems is that parton densities are low and a single parton from the beam 
and single parton from the target participate in a hard reaction. As a result, hard cross 

sections are linear functionals (convolutions) of the appropriate parton densities in the pro- 

II 

jectile and target 'IJ . For instance, once the unintegrated gluon density of the target proton 
is determined from the deep inelastic scattering (DIS) structure function, it would allow a 
consistent description of all other small-x, i.e., high-energy, processes of hard production 
off free nucleons. In contrast to that, in hard production off nuclei the contributions of 
multigluon exchanges with the nucleus are enhanced by a large size of the target. The 
principal consequence is a dramatic breaking of the conventional linear fc^-factorization for 
lard processes in a nuclear environment which, according to the recent extensive studies 
3 fll fl a. be re^aced by a „o„li„ea. ..-fac.o..atio. Namely, one ea. ta.e 

diffractive dijet production p, 1^ as a reference process for the definition of the collective nu- 
clear unintegrated gluon density. Then, it turns out that the so-defined nuclear glue furnishes 
the familiar linear fc^-factorization description of the nuclear structure function F2a{x,Q'^) 
and of the forward single-quark spectrum in DIS (although the linear /c^-factorization prop- 
erty of both observables is rather an exception due to the Abelian feature of the photon). 
Furthermore, the dijet spectra in DIS and single-jet spectra in hadron-nucleus collisions ad- 
mit a description in terms of the same collective nuclear gluon density, albeit in the form 
of highly nonlinear quadratures. The universality classes introduced in js, allow to re- 
late the nonlinearity properties of final states from different partonic pQCD subprocesses 

to the pattern of color flow from the incident parton to final-state dijet. A full derivation 

n 

of nonlinear fc_|_-factorization for all high-energy single-jet spectra was published in [^]; the 
forward quark-antiquark dijet production in DIS and pion-nucleus collisions was studied 
m 3 and 0, respectively; the results for the two-particle pectrum of open heavy flavor 
production g QQ in gluon-nucleus collisions - the dominant source of charm in proton- 
nucleus collisions - were presented in jslOl; quark-gluon dijets in quark-nucleus interactions 
- the dominant source of forward dijets in the proton fragmentation region of proton-nucleus 
collisions - were treated in j^. 

In this communication we report the derivation of nonlinear /c^-factorization for the 
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last missing pQCD subprocess - the production of hard gluon-gluon dijets in gluon-nucleus 
coUisions when the nuclear coherency condition x < xa ~ 0.1 • A^^^^ holds (for the definition 
of xa for a target nucleus of mass number A see below). At the not so high energies of the 
Relativistic Heavy Ion Collider (RHIC), the coherency condition can only be met in the 
proton fragmentation region of pA collisions, where the contribution from gluon-gluon dijets 
is marginal (IS] and references therein). This subprocess will be a principal building block of 
the pQCD description of mid-rapidity dijet production in pA collisions at the Large Hadron 
Collider (LHC), however. The non-Abelian intranuclear evolution of gluon-gluon dijets is 
quite involved - at arbitrary number of colors Nc two gluons couple to seven irreducible 
representations. Based on the reduction of the dijet production problem to interaction with 
the nuclear target of color-singlet multiparton states 0,0, [iQ, U, 3|, we report an explicit 
form of evolution matrices for arbitrary Nc- We demonstrate how the forbidding case of 
seven-channel non-Abelian evolution equations can be diagonalized in an explicit form in 
the large- iVc approximation - this is reminiscent of our finding of the reduction of the three- 
channel non-Abelian evolution for quark-gluon dijets to a two-channel problem ^. 

The production of gluon-gluon dijets on nuclear targets in a limit of strong ordering of 
the rapdities of the produced gluons has been discussed earlier by several authors HQ. 
The starting points are similar, but these works stopped short of the explicit diagonalization 
of their version of the non-Abelian evolution for the four-parton state. Baier et al. reported 
some numerical results for the equal-momentum dijets. As we commented in j^, our ana- 
lytical results ariticipated the trends of nuclear decorrelation of dijets reported in In 



contrast to 



14l |. our results for the gluon-gluon dijets in all color representations are 



presented in the form of explicit quadratures in terms of the collective nuclear unintegrated 
glue and do not assume a stron^_ordering of the gluon rapidities. They fully confirm our 
concept of universality classes 5|, |6[. For instance, the nonlinear k± factorization properties 
of excitation of digluons in higher color representations are the same as those in excitation 
of color-octet quark-antiquark dijets in DIS and quark-gluon dijets in higher color repre- 
sentations in qA collisions. The only difference is in the collective nuclear glue - different 
pQCD subprocesses pick up different components of the color-density matrix for nuclear 
glue. This is also the case for excitation of dijets in the same color representation as the 
incident parton: g {gg}8, 9 {^^js? Q ^ {o'S'js- The diffractive excitation of digluons 
in the antisymmetric octet is similar to diffractive excitation of color-triplet qg dijets in 
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qA collisions and color-octet quark-antiquark dijets in gA collisions. In both g gg and 
q qg processes coherent diffractive excitation of incident partons with net color charge is 
suppressed by a nuclear absorption factor which can be identified with Bjorken's gap survival 
probability |l,5|. 

The further presentation is organized as follows. We start with the discussion of the 
reaction kinematics and the master formula for the dijet cross section in Sec. 2. The 
interaction properties of the two-gluon and three-gluon states are presented in Sec. 3. The 
technically rather involved derivation of the nuclear S-matrix for the four-gluon state is the 
subject of Sec. 4. In Sec. 5 we report the linear fc^-factorization formula for the gluon-gluon 
dijet cross section for the free-nucleon target. The principal new results of our study - 
nonlinear fc^-factorization formulas for gluon-gluon dijets in different color representations, 
their classification in universality classes and a comparison to other dijet processes - are 
reported in Sec. 6. In Conclusions we summarize our main results. 

The technicalities of the construction of the irreducible representations for the two-gluon 
states at an arbitrary number of colors Nc are reported in Appendices A,B. The exact 
integration of non-Abelian evolution equations for the four-gluon system and the derivation 
of explicit quadratures for the dijet spectrum is possible only for large Nf.-, although the 
calculation of higher order terms of 1/Nc perturbation theory is not a problem On the 



and in Appendix 



other hand, the single-jet problem can be solved exactly at arbitrary N^. 
C we show how the coupled seven-channel equations can be exactly diagonalized in the t- 
channel basis appropriate for the single-jet problem. In Appendix D we give a summary of 
different components of the color-density matrix for nuclear glue which enter the description 
of different pQCD subprocesses. 

II. THE MASTER FORMULA FOR GLUON-GLUON DIJET PRODUCTION OFF 
FREE NUCLEONS AND NUCLEI 

A. Kinematics and nuclear coherency 

Our exposition of the master formula for dijet production follows closely our recent work 
on quark-gluon dijets j^. 

To the lowest order in pQCD the underlying subprocess for gluon-gluon dijet production 
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FIG. 1: The rapidity structure of the radiation of gluons by gluons g ^ gg in the nuclear coherency 
region of pA colhsions. 



in the proton fragmentation region of proton-nucleus collisions is a collision of a gluon g* 
from the proton with a gluon qn from the target, 

9*9n 99 ■ 

It is a pQCD Bremsstrahlung off a gluon tagged by the scattered gluon. We do not restrict 
ourselves to the emission of slow, z <^ 1 gluons. In the case of a nuclear target one has to 
deal with multiple gluon exchanges which are enhanced by a large thickness of the target 
nucleus. 

From the laboratory, i.e., the nucleus rest frame, standpoint it can be viewed as an 
excitation of the perturbative \gg) Fock state of the physical projectile \g*) by one-gluon 
exchange with the target nucleon or multiple gluon exchanges with the target nucleus. Here 
the collective nuclear effects develop, and the frozen impact parameter approximation holds, 
if the coherency over the thickness of the nucleus holds for the gg Fock states, i.e., if the 
coherence length is larger than the diameter of the nucleus 2Ra, 



Ir 
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iQ*) 



Mi 
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> 2Ra: 



where 



M2 = ^ + M 



(1) 



(2) 



is the transverse mass squared of the gg state, Pi 2 and zi^2 are the transverse momenta and 
fractions of the the incident gluon's momentum carried by the outgoing gluon one and gluon 



two, respectively {zi + Z2 = 1). The virtuality of the gluon g* equals (Q*)^ = (p*)^, where 
p* is the transverse momentum of g* in the incident proton (Fig. Q). In the antilaboratory 
(Breit) frame, partons with a momentum xp^ have a longitudinal localization of the order 
of their Compton wavelength A = l/xp^, where p^ is the momentum per nucleon. The 
coherency over the thickness of the nucleus in the target rest frame is equivalent to the 
spatial overlap of parton fields of different nucleons at the same impact parameter in the 
Lorentz-contracted ultrarelativistic nucleus. In the overlap regime one would think of the 
fusion of partons form different nucleons and collective nuclear parton densities Q]. The 
overlap takes place if A exceeds the Lorentz-contracted thickness of the ultrarelativistic 
nucleus, 

X = ^>2R^.^, (3) 

xpN Pn 

which is identical to the condition (jT]). 

Qualitatively, both descriptions of collective nuclear effects are equivalent to each other. 
Quantitatively, the laboratory frame approach takes advantage of the well developed 
multiple-scattering theory of interactions of color dipoles with nuclei . From 

the practical point of view, the coherency condition x < Xa restricts collective effects in hard 
processes at RHIC to the proton fragmentation region of pA, dA collisions, but at LHC our 
treatment will hold down to the mid-rapidity region of pA collisions. The target frame rapid- 
ity structure of the considered g* gg excitation is shown in Fig.[T] The (pseudo)rapidities 
of the final state partons must satisfy 771,2 > 'Ua = log l/x^- The rapidity separation of the 
two hard gluon jets, 

= log ? ' (4) 

is considered to be finite. Both jets are supposed to be separated by a large rapidity from 
other jets at mid-rapidity or in the target nucleus hemisphere; the gaps between all jets, 
beam spectators and target debris are filled by soft hadrons from an underlying event. 

Clearly, the incident gluon g* is a parton of the color-singlet beam hadron and is accom- 
panied by comoving spectators. However, the effect of interactions of comoving spectator 
partons cancels out upon the fully inclusive integration over the spectator phase space [ill , 
an explicit demonstration of such a cancellation is found in jj]. The properties of the beam 
hadrons only define the longitudinal and transverse momentum spectrum of g* in the beam, 
and the problem we treat here will be a building block of the pQCD description of mid- 
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rapidity to proton hemisphere dijets in pA colhsions at LHC. 
B. Master formula for excitation of gluon-gluon dijets 

In the nucleus rest frame, relativistic partons g*,gi and g2, propagate along straight-line, 
fixed-impact-parameter, trajectories. To the lowest order in pQCD the Fock state expansion 
for the physical state \g*)phys reads 

\g*)phys = \g*)o + '^{zi,r)\gg)o, (5) 

where \Ef(zi,r') is the probability amplitude to find the gg system with the separation r 
in the two-dimensional impact parameter space, the subscript "0" refers to bare partons. 
The perturbative coupling of the g* gg transition is reabsorbed into the lightcone wave 
function \E'(2;i,r). We also omitted a wave function renormalization factor, which is of no 
relevance for the inelastic excitation to the perturbative order discussed here. The explicit 
expression for \E'(zi,r) in terms of the gluon-splitting function and the gluon virtuality 
(Q*)^ will be presented below. For the sake of simplicity we take the collision axis along the 
momentum of the incident quark g*, the transformation between the transverse momenta in 
the (7*-target and p-target reference frames is trivial jj] . 
If b is the impact parameter of the projectile g*, then 

bi = b - Z2r, b2 = b + zir . (6) 

By the conservation of impact parameters, the action of the S-matrix on \a)phys takes a 
simple form 

S|/W = Sg{b)\g*)o + Sg{bi)Sg{b2)^{zur)\gg)o 

= Sg{b)\g*)phys + [Sg{b^)Sg{b2) " S g (b)]^ {z, , \ g g) ^ - (7) 

In the last line we explicitly decomposed the final state into the (quasi)elastically scattered 
\g*)phys and the excited state \gg)o. The two terms in the latter describe a scattering on the 
target of the gg system formed way in front of the target and the transition g* — > gg after the 
interaction of the state |5'*)o with the target, as illustrated in Fig. |21 The contribution from 
transitions g* gg inside the target nucleus vanishes in the high-energy limit oi x <xa^- 

^ In terms of the lightcone approach to the QCD Landau-Pomeranchuk-Migdal effect, this corresponds to 
the thin-target hmit (12j . 
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a) b) c) 

FIG. 2: Typical contribution to the excitation amplitude for gA — > gig2X, with multiple color 
excitations of the nucleus. The amplitude receives contributions from processes that involve inter- 
actions with the nucleus after and before the virtual decay which interfere destructively. 



We recall, that the s-channel helicity of all gluons is conserved. 

The probability amplitude for the two-jet spectrum is given by the Fourier transform 

I d^b,£b2exp[-^ip,b^+p,b2)][Sgib^)Sgib,) - S,(6)]vl/(^i, r) (8) 

The differential cross section is proportional to the modulus squared of (jH)), 

/ d%d%exp[z{p,b[+p,b'MSl{K)Sl{y,) - Slib'Wiz^r') 

X I d^bid^b2exp[-tip,bi+p,b2)][Sg{bi)Sg{b2)-Sgib)]^{zi,r). (9) 



The crucial point 



antiparton 



IG 



is that the hermitian conjugate can be viewed as the S-matrix for an 



Consequently, the four terms in the product 
[S,{b[)S,ib',) - S,ib')nS,{b^)S,ib2) - S,(6)] 



admit a simple interpretation: 



S%,{b',b) = St(6')S,(b) 



(10) 



can be viewed as an S-matrix for elastic scattering on a target of the g*'g* state in which the 
(anti)gluon g*' propagates at the impact parameter 6'. The averaging over the color states 
of the beam parton g* amounts to the dipole g*g* being in the color singlet state. Similarly, 

:(3) 



S'Xjb',bub2) = St(b')S,(6i)S,(b2) 



S%Jb,b[,b',) = St(b;)St(6',)S,(6) 



9l929 

^Z'.9.9.(^'^^^2MM) = sl{b[)Sl{b',)SMS,{b2) . 
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(11) 
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q(4) ^ ^(2) 

^9i929^92 ^99' 






q(3) 
^99i92 



9 9i92 



FIG. 3: The S-matrix structure of the two-body density matrix for excitation g 9192- 



describe elastic scattering on a target of the overall color-singlet three- and four-gluon states, 
respectively. This is shown schematically in Fig. El Here we suppressed the matrix elements 
of S^"^ over the target nucleon, for details of the derivation based on the closure relation, see 
P]. Specifically, in the calculation of the inclusive cross sections one averages over the color 
states of the beam gluon g, sums over color states X of final state gluons qi,g2, takes the 
matrix products of and S with respect to the relevant color indices entering S*-""-* and sums 
over all nuclear final states applying the closure relation. The technicalities of the derivation 
of S^"-* will be presented below, here we cite the master formula for the dijet cross section, 
which is the Fourier transform of the two-body density matrix: 

da{g* gig2) 1 



X exp[-ip2(&2 - fa'a) - iPiibi - 61 - b2W{zi, b[ - b'^) 

E(^l{sSk....(^i' b,) + Sf4.(b', b) - S(f^,^,(6, b[, b',) - Sf4,,(6', 6„ b.)} 



(12) 

in) 



Hereafter, we describe the final state dijet in terms of the jet momentum p = p^, z = zi, 
and the decorrelation (acoplanarity) momentum A. = p^ + P2. We also introduce 



s = 62 - 62 , 



(13) 



in terms of which bi — b[ = s + r — r' and 

exp[— 2^2(^2 — ^2) ~ Wi{bi — b\)] = exp[— ?As — ipr + ipr'] 



(14) 



so that the dipole parameter s is conjugate to the acoplanarity momentum A. 
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III. CALCULATION OF THE 2-PARTON AND 3-PARTON S-MATRICES 



A. The gluon-nucleon S-matrix and the /c^-factorization representations for the 
color dipole cross section 



In order to set up the formalism, we start with the S-matrix representation for the cross 
section of interaction of the gg color dipole with the free-nucleon target. To the two-gluon 
exchange approximation, the S-matrix of the gluon-nucleon interaction equals 

S^(6) = 1 + iT'^VaXib) - ]^T^T-x\b) , (15) 

where T°- is the SU{Nc) generator in the adjoint representation, {g^\T"-\g^) = —ifabc, and 
T°-Vax{b) is the gluon-nucleon eikonal for single gluon exchange. The vertex Va for excitation 
of the nucleon g^-N N* into a color octet state is so normalized that after application of 
closure over the final state excitations A^* the vertex g°'g^NN equals (A^| V"j^Vf,|A^) = 5ab- The 
second order term in p5|) already uses this normalization. The S-matrix of the gg-rmcleoii 
interaction equals 

.(2).. .^ (iV|Tr[S^(bi)SUb2)]|iV) 

\, (bu &2) ^N\Trl\N) " ^^^^ 

The corresponding profile function is r2(6i, 62) = 1 — ^fgi^i, &2)- The dipole cross section 
for interaction of the color-singlet gg dipole r = bi — 62 with the free nucleon is obtained 
upon the integration over the overall impact parameter 

a{r) = 2 1 £b,T2ib,, b,-r) = CAj d^b, [x(&i) - x(&i - r)]\ (17) 

where Ca = T°'T°' = is the gluon Casimir operator. Eq. [T7I sums up the contributions 
from the four Feynman diagrams of Fig . |3] and is related to the gluon density in the target 
by the /c^-factorization formula [lol bcj l 

(r{x,r) = -^J d'^K,f{x,K)[l -ex.p{iKr)], (18) 

where Cp = {N^ — l)/2Nc is the quark Casimir. Recall that the unintegrated gluon density 

a log 
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FIG. 4: The four Feynman diagrams for the gluon-gluon dipole-nucleon interaction by the two- 
gluon pomeron exchange in the t-channel. 



was defined with respect to the qq color dipole probe, and is related to /(x, k) through 

/(x, k) = ■ K ) . (20) 



Hereafter, unless it may cause confusion, we suppress the variable x in the gluon densities 
and dipole cross sections. The energy dependence of the dipole cross section is governed by 
the color-dipole leading Log^ evolution 3, 121I , the related evolution for the unintegrated 
gluon density is described by the familiar momentum-space BFKL (Balitsky-Fadin-Kuraev- 
Lipatov) equation j^]. 

The S-matrix for coherent interaction of the color dipole with the nuclear target is given 
by the Glauber-Gribov formula 0, 



where 



S[6,a(r)] = eM-\<r)T{b)] 



/oo 
-00 



(21) 



(22) 



is the optical thickness of the nucleus. The nuclear density n /i (&, rS] is normalized according 



In the specific case of Sf}„*{b' , b) the color dipole equals 



99 



b — b' = s + zr — zr' 



(23) 



and S^gJg*{b' , b) entering Eq. (fT!?jl will be given by the Glauber-Gribov formula 



S%{b',b) = S[b,ais + zr-zr')]. 



(24) 
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3 

a) 



b) 



FIG. 5: The color dipole structure of (a) the generic 3-gluon system of dipoles and (b) of the g'gig2 
system which emerges in the S-matrix structure of the two-body density matrix for excitation 
9 9i92- 

B. The S-matrix for the color-singlet ggg system 

Here, there are two possibilities to couple three gluons to a color singlet, but only the 
/-coupling is relevant to our problem, see also Appendix C. 

For the generic 3-gluon state shown in Fig. El the color-dipole cross section equals 

62, bs) = \Hr12) + (7(r23) + (7(r3i)], (25) 

where = — hi. The configuration of color dipoles for the case of our interest is shown 
in Fig. El (see the related derivation in For the g'qig2 state the relevant dipole sizes in 

equal 

rg,g*' = bi-b' = s- zr, 
rgi92 = b2-bi = r, 

^5*'52 = b' -b2 = s + r - zr', (26) 

so that 

^9*'gi92 = + a{s + r- zr') + a{s - zr')] , 

^9'9W, = li^i-r') + a{s -r' + zr) + a{s + zr)] . (27) 

The overall color-singlet 3-gluon state has a unique color structure and its elastic scattering 
on a nucleus is a single- channel problem. Consequently, the nuclear S-matrix is given by the 
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single-channel Glauber-Gribov formula 



(28) 



IV. DIPOLE CROSS SECTION OPERATOR FOR FOUR-GLUON STATES 

We now come to the major technical novelty of this paper, the dipole cross section matrix 
for the four-gluon system. Here, as in the previous applications for qqqq and qqgg systems, 
the large Nc limit offers a particularly useful expansion. We thus discuss the general case 
in which the gluon is transforming in the adjoint representation of SU{Nc). In order to 
construct the relevant four-gluon states, our first task is to decompose the product-states 
of the adjoint x adjoint system into irreducible representations. The adjoint (or regular) 
representation of SU (Nc) has A^^ — 1 states, and the Clebsch-Gordan series for the product 
of two adjoint s reads 



In the last line, we named the representations by their SU (3) dimensions, except for one of 
the symmetric representations that vanishes for Nc = 3, and will be referred to as Rj. 

A. Projection operators onto irreducible representations 

The derivation of projectors onto the representations ()29|) is a lengthy, though standard 
I I 

exercise [25]. While in our earlier solution of the analogous problem for the qggq-sjstem 
it had proven convenient to represent gluons in a double line notation as pointlike quark- 
antiquark systems, here we find it expedient to stick to purely adjoint-index tensors. We 
present here a sketch of the construction of irreducible representations (j^^ . the details are 
given in Appendices A and B 



(iV^ - 1) X (iV^ _ 1) 



1 + (iV^ - 1)^ + (iV,2 - l)s 
, (iVc'-4)(iV2-l) , dN! - A){N^ - Ih* 




iV2(iV, + 3)(iV,-l) , iV2(iV,-3)(iV, + l) 



4 4 
1 + 8a + 8s + 10 + To + 27 + i?7 • 



(29) 
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If t*^, a = 1 . . . A*"^ — 1 are 5'f/(A''c)-generators in the fundamental representation, the 
famihar /- and ci-tensors are defined through 

ft' = ^Sab 1 + ^ {dabc + tfabc) f , (30) 

and T^" = —ifabc are the SU{Nc) generators in the adjoint representation. 

First we decompose the product representation space into its symmetric and antisym- 
metric parts, respectively: 

Kd = ^acSbd = S^d + -^cd^ (31) 

where 

S^d = \ {^aAd + SadSbc) , ^ed = ^ {^aAd ' SadSbc) ■ (32) 

The complex 

"^^cd = -^{dadkfkbc + fadkdkbcj (33) 

as well as 

[-^t]cd = dackdkbd , [Dv\^ct = dadkdkbc , [Ds]'^ = dabkdkcd (34) 

also prove helpful. All the above defined tensors S, A, iY, Ds,Dt, Du are hermitian. 

The 5'f/(A^c)-pi'ojectors into the singlet as well as the two adjoint multiplets have mani- 
festly the same form as their well-known = 3 counterparts: 

Pl^td = -j;^^abSed (35) 
-P[8A]cd = J^fabkfkcd = J^^ifabkifkdc (36) 
P[^sTcd = 1 ^ dabkdkcd = ^2^4 i^'^) 

For the construction of higher multiplets a useful quantity is Q"^ = 4 ■ Tr f^tH'f^ , which, 
with indices suppressed, equals 

Q = ^(25 - (AT^ - 1)P[1]) + ^(A + Du - D,) + (38) 



Then, the crucial observation is 



that 



/V2 _ 1 - 4 
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TABLE I: Properties of Multiplets 



symmetric 



antisymmetric 



Name of rep. 1 85 27 R7 8a 10 10 

Dimension 1 - 1 ^.^(^c+3)(iVc-i) iV.^(iV.-3)(iVe+i) ^2 _ ^ iN?-mN?-i) jN^-mN?-!) 

Casimir C2[i?] Nc 2(iVc + 1) 2(iVc - 1) Nc 2Nc 2Nc 

Apparently in the subspaces that are projected onto by 

S^^S- P[l] - P[%s\ , A^^A- P[%a] , (40) 

the operator t — Q'^ — {1 + Q){t — Q) vanishes and Q has eigenvalues ±1. We can thus 
write down projection operators that decompose »S_l,^_l further, as 

Pi = ^ (1 ± Pt = ^ (1 ± Q)^± ■ (41) 

Checking how many states are contained in P^ , , one may identify, that P^ — 
P[10],PJ^ = P[10],P^^ = P[27],PJ^ = P[R7]. In convenient form, with all indices shown 
again, they read: 

^[TO]S = - P[8a]:^ - iY:i) . (42) 



lab 

\cd 



P[^7r^ = ^((iVe + 2)<S:,^-(Ar, + 2)(Ar,-l)P[l]^ 

- i(Ar, - 2)(Ar, + 4)P[85]:^ + Y^mfa + [Du]:^)) (43) 
Pmtd = ^((^^c-2)5,t + (iV,-2)(iV,+ l)P[l]:^ 

+ 1{N. + 2m - m^s]^^ - y([A]»^ + [Du]:^)) . (44) 

It is now a simple matter to obtain the quadratic Casimirs (i.e. color charge squared) of 
the individual multiplets, which can be found in Table I. 

B. Multigluon states 

The above given projectors can be used to construct the color-space wave function of 
the multigluon states relevant for us. The four gluons have to be in a total color singlet, 
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and all possible states are exhausted by coupling a chosen pair of gluons to all possible 
multiplets and the remaining two to an anti-multiplet. The choice of pairs is of course 
arbitrary. Because averaging over colors of the incoming gluon amounts to the initial state 

|in) = , ^ \8a8a), (45) 



a convenient choice is the s-channel pairing 

\RR) ^ \{[g-{b,)^g\b,)]^^[g^{b[)®g\b',)]_}^) 

^ : P[R]t, \g%b,) ® g\b2) ® g%b[) ® g\b',)) . (46) 



/dim[i?] 

The basis of color-singlet four-gluon states which contribute to the 
non-Abelian evolution of gluon-gluon dijets in our problem consists of 
|11), |8a8a), I8585), |10T0), |T010), |2727), {RrRi). The mixed-symmetry color-singlet 
states like 18^85) are possible but decouple from the above states. 

We note in passing, that the single-jet spectrum derives from the dijet spectrum by 
integration over all A, which entails s = 0. For studying the transition from the dijet 
to single-jet problem an alternative, t-channel, pairing of gluons, {gig'i) and {g2g2), proves 
to be a more convenient one. Evidently, the multiparton S-matrices for different choices 
are related by a trivial permutation of the gluon impact parameters. For the reference 
purposes, in Appendix C we demonstrate how the coupled seven-channel problem is exactly 
diagonalized for the color-dipole configuration appropriate to the single-jet problem. 

C. Multiparton S-matrix and the four-body color dipole cross section operator 
for the free-nucleon target 

The frozen impact parameter approximation leads to a four-gluon S-matrix of the form 

sS^)(bi, 62, b[, b',) = S;v(&i) ® 5^(62) ® S^(6;) ® S^ib',) . (47) 

Then, upon using eq. (fT5|l . on a color-singlet four-particle state, the S-matrix S^'' takes the 
form 

Sg)(bi, 62, b[, b',) = 1 - ^[Tfx(bi) + T»x(b2) + rf,x(b;) + TMb2)]' 
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-T^T^xipi)x{b2) - n>nxib'Mb'2) 
-T^T^,x{b2)x{b\) - T^T^,x{bMb'2) 

-T,^T^,x{bi)x{b[) - T^T^,x{b2)x{b'^) ■ (48) 

Here products like Tf Tg" are shorthands for 1 <8) 1 <8) T"", what acts in the space spanned 
by gluon states gl® g^® gy ® g^ ■ 

In the following we shall heavily exploit the color-singlet condition for the four gluon 
system, which reads 

(Tf + + Tf, + T2", )\RR) = Q. (49) 

As we work with states in which the dipole 12 and the conjugate dipole 1'2' are in definite 
color representations, we have 

(Tf + T^f\RR) = {T^, + T^,f\RR) = C2[R]\RR) , (50) 

which we can use to simplify the cross product terms, e.g. 

T^T^ = T^,T^, = ^{C2[R\ - 2Ca) , (51) 

where we used that T^T^ = Ca for all j ^. Now notice that for the operators 

Tg = Tf + = T" ® 1 ® 1 ® 1 + 1 ® T" ® 1 ® 1 
Tg, = Tf, + T2'5 = 1 ® 1 ® ® 1 + 1 ® 1 (8 1 T'^ 
= - T2" = T" 1 ® 1 1 - 1 ® T" 1 ® 1 
A^, = Ti'^ - T2", = 1 1 (g) T" (g) 1 - 1 (g) 1 (g) 1 T" (52) 

we have 

T^Al = (T^)^ - {T^f = = rg,A^, , (53) 
and, because of the color singlet condition = — T^/, also 

rS,A^ = T^A^, = 0, (54) 



Here a sum over a, but not j is implied 

17 



and, effectively 



rpa rpa 



rpa rpa 



(55) 



Now insert 
1 



= -{n + Ad) ; ^2 = \iT^D- A^) ; T^, = ^(T^, + A^,) ; T,^, = 1(T^, - A%,) , (56) 
into PHjl . and use relations (jSH) and (jSEI) to obtain 

1 - s^^\b„b„b[,b',) = ^CA{xHb,) + xHb,) + xHb[) + xHK)) 

+^{C2[R]-2CA){x{bi)x{b2) + x{b[)x{b'2)) 
-^C2[R]{x{b[)x{b2) + X{bi)x{b'2) + x{bi)x{b[) + x{b2)x{b'. 



4 
1 

i 

4 



We can now go ahead and complete the squares, to obtain 



A%A%(xibMb[) + Xib2)xib',) - Xib[)xib2) - x{b,)xib'. 



(57) 



l-S^^\b^MAA) = -^{2Ca-C2[R]) 



+-C2[R] 



(x(6i) - X(&2)) + {x{b\) - x{b',) 
{x{b,) - xib\)f + {xib,) ~ x{b'. 



+ (x{b\) - xib2)f + {xib,) - xib',. 



J A« A" 



X{b[)-X{b2)) + (x{bi) - X{b'2 



-(X(bi)-x(b'i))'-(x(b2)-x(&2) 



(58) 



Finally, using ()17j) we obtain the following form for the dipole cross section operator for 
the four-gluon system 

{R'R\a^^\b,, 62, b'l, b',)\RR) = A«5k',r ■ [a{b^ - b^) + a{b\ - b',] 
, (1 - Aij) 



-Sr' 



R'R ■ 



a{bi - b[) + a{b2 - 62) + a{b\ - 62) + (y{bi - b'2 



Here the parameter 



1 - 



cm 

2Ca 



(59) 



(60) 
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FIG. 6: The matrix element {RjRj\0\RiRi) . 



enters the diagonal part of the cross section and can be found for individual multiplets in 
table I. In the off-diagonal piece we introduced the combination 

62, K, b'^) = a{b[ - 62) + (T{b, - b'^) - or(6i - b\) - a{b2 - b'^) 

= a{s + r) + a{s — r') — a{s) — a{s + r — r') , (61) 

the same structure of dipole cross sections made already an appearance in our previous 
solutions of the qqqq and qgqg dipole cross section matrices 0,0]. Eq.(|^ is the central result 
of this subsection. We now turn to the evaluation of the matrix elements {R' R'\A'}jA'^,\RR) . 



D. Evaluation of the off diagonal matrix elements 

For R ^ R', we have 

{R'W\a^'\b„ 62, b[, b',)\RR) = JR'R'I^D^URR) ^^^^^ y^^ (g2) 

We recall, that 

A^A^, = (t'^ ® 1 - 1 (g) T") O (t'^ ® 1 - 1 O T") 

= ® 1 ® T'* ® 1 - 1 ® ® T'' ® 1 - ® 1 ® 1 ® + 1 ® 1 ® . (63) 

One can easily convince oneself that this is really a transition operator, i.e., it has only matrix 
elements between different multiplets, and even more they must be of different permutation 
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TABLE II: Matrix elements of the Operator O 



A \N,{Nl - 1)2 
P[U lN,{Nl - 1) 
iY 



^[1] ^^[85] 

N, |Ar,(iV2 - 1) 

N, \N,{N^, - 1) 





i(iV2-l)(Ar2-4) 





-l(iV2-l)(iV2-4) 
i(iV2-l)(Ar2-4) 




symmetry, i.e. the transitions are between symmetric and antisymmetric multiplets. One 
can then use permutation symmetry, to obtain, effectively, 

A^A" , = 4 ■ (1 ® T" ® 1 ® T") 

= -4 • (1 ® T" ® 1 O (T")*) = 4 • 6 , (64) 

which means 

{R'R\a^^\bi, 62, K, b'^)\RR) = --^ ■ {R'R\d\RR) ■ 62, K, b'^) . (65) 

Between four gluon states with color wavefunctions v4^", BJ^^ the matrix element of the 
operator O is evaluated explicitly as 

{A\0\B) = ^fnan'^fall' ^^^'^nn'^C , (66) 

see also Fig.®. 

In table II we collect the contractions of the operator O between a convenient choice 
of tensors, from which the matrix elements of the relevant four gluon states (jlUj) can be 
reconstructed. The tensor iY decouples completely - all its off-diagonal elements vanish. 



E. The dipole cross section operator for the four gluon system 

We now come to the final result for the four-body dipole cross section operator in the 
s-channel basis described in Sec. IV B above, in which the two gluons from the amplitude 
and from the complex conjugate amplitude, respectively, are in definite color multiplets. 



1. Diagonal elements 

The diagonal elements for the four-gluon system are even simpler than for the quark- 
antiquark-gluon-gluon system studied in and are expressed in terms of two combinations 
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of color-dipole cross sections: 



El = a{b, - 62) + ^{b\ - b',) = a{r) + a{-r') , 
T = a{b, - b[) + (7(62 - b'^) + a{b[ - 62) + - b'^) 

= a{s + r) + a{s — r') + cr(s) + cr(s + r — r') . (67) 

Making use of table I, we find 

(8^8^|a(^)|8A8^) = (8585|^(')|8585) = ^r + isi 

(10lO|a(^)|10TO) = (I010|a(^)|I010) = \r = ■ \r 

(4) 2(iV, + l) 1 1 C2[27] 1 1 

(27271^^12727) = " 4- " = ^ " 4- - 

(i?,i?,|a \RrRr) - ■ -r + = ^— ■ -r + (68) 

We recall that in the limit of r = r' = the gluon pairs collapse into pointlike partons in the 
color representation Ri. In this limit Si = 0, and the diagonal matrix elements are simply 
proportional to the Casimir operators as it must be 0|: 

(i?.i?.|a(^)|i?.i?.) = %t^a(5). (69) 

Notice that the matrix elements (2727|a(^) |2727) and {R7Rj\a^^'^\R7Rj) are related by the 
transformation A^^^ — > —Nf.-, for the discussion of a similar symmetry in the quark-gluon dijet 
production see Ref. j6| 

2. Off- diagonal elements 

Making use of Tables 1, 1 1 we have 



(ll|a(^)|8^8^) = ^^.(7(bi,b2,6;,6'2) 

In^ - 1 



(85851^(^)18^8^) = -^■r](6i,62,6'i,&2) 



(2727|aW|8^8^) = ■ l](bi, 62, b',, 6^) 
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{RjRr\a^'^\8A8A) = -^\/ ' ^(^i' ^2, b'l, ^2) 
(858s|a(^)|10lO) = ■ fi(6i, 62, K, b',) 

= (8585|a^^^|T010) 



(4)1 



(2727|<T^^^|1010) 



(70) 

We again observe the curious symmetry the matrix elements (2727|(j*^'^'' |8/i8a) and 
{RtRyIct^'^^ISaSa) are related by the transformation — * — A^c, the same is true of the 
matrix elements (2727|<tW|T010) and (i?7i?7|(j(^)|10T0). 



F. Large-A^c properties of the dipole cross section matrix 

In conjunction with the Glauber-Gribov form for the nuclear S'-matrix, the dipole cross 
section operator a^^^ solves the problem of non-Abelian intranuclear evolution of the four- 
gluon system. Being a symmetric matrix, the dipole cross section operator could readily be 
brought in diagonal form, and the Fourier-transform could finally be performed numerically. 
In practice, however one would encounter a number of obstacles when proceeding along 
these lines. First, for the case of the two-particle inclusive spectrum, the eigenvalues of 
the dipole cross section operator will be non-algebraic functionals of the free-nucleon cross 
section operator, and second the dipole cross section itself has a non-analytic dependence on 
dipole size, which would ultimately determine the asymptotics of the Fourier-transforms. It 
is therefore convenient that the large- iVc expansion offers a path to analytic formulas, which 
can be interpreted in a transparent way. 

We start from the observation, that at large- A^c the cross section operator a^'^^ assumes a 
block diagonal form. Apparently, transitions between representations which have dimensions 
that grow with the same power of Nc are parametrically of order N^, whereas transitions to 
the next-larger /smaller block are suppressed by A""^. 
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In the space of four parton-states we introduce the projectors 

Vi =|11)(11| 

= |8a8a)(8a8a| + 1 8585) (8585 1 

= |10T0)(10T0| + |T010)(T010| + |2727)(2727| + \R7R7){R7R-j\ , (71) 

which aUow us to isolate the blocks Vi a^'^^ Vi, the first one being a one- by-one matrix: 

Via^^^Vi = El . (72) 

The vector |ei) = |11) can be viewed as an eigenvector of block 1 with eigen-cross-section 
Si. 

The second block is written, in the two-dimensional subspace of octets, as 

i[T + 2Ei] -iQ ^ 

-iQ i[T + 2Si]^ ■ 

Its eigenvectors are 



(73) 



|e2) =^(|8a8a) + |8585)) 



(74) 



and belong to the eigen-cross-sections 



El + -[r - 1]] = - (El + a(6i - b'^) + (7(62 - b'^ 



2^ 



- (7(r) -I- (j(— r') -I- (t{s) -\- a{s + r — r') 
^Ei + ^[r + = ^ [Ei + a(6; - 62) + (t(6i - 6^; 



4' ■ 2 

= ^ o-(r) (j(-r') + (7(s r) (t(s - r') . (75) 

Finally, the third block accounts for the multiplcts that have 0{N^) states. Here we no- 
tice that at large all higher multiplets interact as two color-uncorrelated gluons. For 
instance, the Casimir operators for these multiplets approach C2{Ri) — 2Ca and the di- 
agonal cross-sections become identical. In matrix form, where the rows refer to states 
1 1010), |T010), 1 2727), \R7R7) we write 







-\n -\n 



(76) 
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Its eigenvectors are easily seen to be 



|e4) = ^(|1010) + 11010) + 12727) + \R^R^ 
leg) = J(|10T0) + |T010) - |2727) - |i?7i?7 



2 



ke) 



i=(|1010)- 11010) 

1 

71 



|e7) = 4^(12727) -|i?7i?7)), (77) 



with eigenvalues 



\{r-n) = a{h-b[) + aib2-b',) 



= (T[s) + (T[s + r — r ) 

Ss = ^(r + n) = a{b\ - &2) + a(6i - b',) 

= a(s + r) + a(s — r') 

S6 = S7 = ir = ^(S4 + S5). (78) 

We finally observe, that the eigenstate jee) decouples exactly from our problem, which 
is seen readily from the summary of the off-diagonal elements fj70p . Couplings between the 
above diagonalized matrix blocks are of 0{N~^). In the basis of eigenstates |ei), . . . , ley) we 
can collect the A^c-suppressed off-diagonal elements as 

uj{s,r,r') = --^j^n{s,r,r')[\ei){e2\ + |ei)(e3| + |e4)(e2| - |e5)(e3| +h.c.} . (79) 

It is easy to check that the off-diagonal elements containing the state ley) are 0(A''~^) 
and this state decouples at 0{N^^). At 0{N^^) there are also corrections to the diagonal 
matrix elements, in the sector of large symmetric representations 27, Rj. They are however 
not relevant for our problem, the perturbation theory treatment of such corrections is 
found in |2]. The crucial observation is, that in the inclusive dijet spectrum, with summation 
over all colors of final state gluons, the final state projection simplifies to 



^(X| =5]v^dim[i?](i?/?| =^ + iV£y2(e^+ . (80) 

^ ^ 1 8a+8s 10+T0+27+R7 

If we remember, that averaging over incoming colors shall introduce another factor 1/Nc, 

|in) = ^=i==|8^8^) = ^(|e2) + ks)) , (81) 
/dim [8] V^N^ ^ ' 
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we see that the large number of states in higher multiplets can overcome the suppression 
of their excitation. Furthermore, the excitation of singlet states is large- A^c suppressed and 
will be neglected from now on. 

To develop the large- A^^^ perturbation theory, we decompose the free-nucleon cross section 
operator as 

a(^) = S(°)+c^, (82) 

where 

S(o) = + P2a(^)P2 + r^^a^^'^V^ = S,|e,)(e,| (83) 

i=i 

is the block matrix that is diagonalized by the basis |ei) . . . ley). Now, the nuclear S'-matrix 
is obtained from the free nucleon dipole cross section by means of the Glauber-Gribov 
exponentiation 

S[b,a(^)(s,r,r')] =exp[-la(^)(3,r,r')T(6)]. (84) 

To the first order in the off-diagonal perturbation ci) we can easily establish 

S[6, S[6, = --T(6) /' d(3S[h, (1 - (3)t^^^] CoS[h, /5E(°)] + 0{N-^). (85) 

2 JO 

V. LINEAR fc^-FACTORIZATION FOR DUETS FROM THE FREE NUCLEON 
TARGET 

At this point we are in a position to give our result for the process g*gN 9x92 on the 
free nucleon target. After integrating over the overall impact parameter, our master- formula 
assumes the form 

xY,{X\t(s,T,r-)\m) (86) 

X 

We evaluate the cross section for an incoming gluon, which entails an average over its 
incoming colors and the initial state \in) of Eq. ()45|) . Likewise in the final state we sum over 
all color states. Then, the calculation of the inclusive cross section involves the evaluation 



25 



of the following matrix elements of the four parton cross section operator 



E 



, ^ dim[i?] 
1 1 

XI - +— T + 



{RR\a^^\s,r,r')\8A8A) = ^{s,r,r') 



N^{N, + 3){N,-1) 1 N, + 3 
^' ' \ 4(iV2 _ 1) 2iVe V iV, + 1 
8s ""^ T ' ' ' 

^1 27 



+ 



iV2(iV,-3)(iVe + l) 1 N,-3 ' 
\ 4(iV2-l) 2iV,ViV,-l^ 
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iJ7 



= if](s,r,r')(4^+ ^ ). (87) 

8s l+27+iJ7 

Here we indicated the contributions from excitation of individual multiplets. Q{s, r, r') is 
the same as in eq. ()61|) . now in the relevant coordinates, explicitly 



r2(s, r, r') = J d^Kf{K) exp[iKs](l — exp[iK;r])(l — exp[— i/tr']) . 



(88) 



The diagonal contribution from final state gluons in the antisymmetric octet is easily con- 
structed from the results of sec. IV: 

(8a8a|S(s, r, 018^8^) = a(3)(6i, 62, h') + a^^\b\, b'„ b) - a(s + ^(r - r')) 
-{8ASA\cr^^\s,r,r')\8A8A) 



= - cr(s — zr' + r) + a{s — r' + zr) — a{s + z{r — r')) — a{s + r — r') 

+cr(s — zr') + cr(s + zr) — a(s + z{r — r')) — a{s) 
1 

~4 

d?Kf{K) exp[iKs] 
x< -(exp[iK;r] — exp[i2;K;r] ) ( exp[— i/tr'] — exp[— i^;^^' 



(j(s — r') + (t(s + r) — a(s + r — r') — a(s) 



■!)(' 

+ -(^1 — exp[i2K;r]^ (^1 — exp[— 22;/tr' 
-(^1 — exp[iK;r]^ ~ exp[— ?Kr']^ | . 



(89) 
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The dijet cross section for production of gluons in the antisymmetric octet reads therefore 
daNig* {gi92}sA) 



^/(A) { \^{z,p-A)-^iz,p-zA)\' + 



dzd?Ad?p 2(27r)^ 



]-\m{z,p)-m{z,p-A)\^}, (90) 



2' 

The wave function of the gluon-gluon Fock state of the physical gluon enters our analysis as 
the recurrent quantity 

- *(^.p - = 2P„M - (p,",)."^,. )'. (91) 

where Pgg{z) is the familiar gluon splitting function, 



1 — z z ,^ , 

+ z + z(l-z) 



z 1 — z 

= z{l-z){Q*Y. (92) 

Now notice, that the last contribution in (j90|) would be canceled exactly, in the fully 
inclusive sum, by the excitation of symmetric octets. Summing over all possible final states 
we end up with 



+ \^{z,p) - ^!{z,p- z\)\^ 
+ \^{z,p)-^>{z,p-A)\^}, (93) 



where now the last line sums up the contribution from excitation of 1, 27, -R7, whereas the 
first two terms represent the sum of octet, 8^,85 final states. No large approximation 
has been invoked here, what would change with Nc is only the composition of the final state, 
where the excitation of color singlet states is 0{N^'^). The absence of decuplet excitation 
is a consequence of the excitation mechanism being single-gluon exchange. Our result eq. 



further elucidating 



(inni) is of course nothing but the differential version of eq. (82) of Ref. 
the color-composition of the final state. By itself it would find interesting applications to 
dijet production in a regime where saturation/absorption effects are not strong, e.g. central 
dijets at HERA or Tevatron. 
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VI. NONLINEAR fc_L-FACTORIZATION FOR DUETS FROM NUCLEAR TAR- 
GETS 



A. Dijets in color-octet final states 

In order to evaluate our master formula for the case of octet-final states we have to 
collect the various multiparton ^-matrices. The three- and two-body ^-matrices are single 
channel problems and can be written in terms of the Glauber-Gribov exponential as 

sS^r = ^^W.A = S[fe, \<r)] S[6, \a{s - zr' + r)] S[6, \a{s - zr')] 

-(3) 



S^l^^,. = S[6, = S[b, ia(r')] S[6, - r' + zr)] S[b, + zr)] 



S%, = S[b, a,*,,,] = S[b, ais + z{r - r'))] = S'[b, ^ais + ^(r - r'))] (94) 



Notice that the factorization of ^-matrices is an exact consequence of the form of the dipole 
cross section for the three-gluon state and of the Glauber-Gribov exponentiation valid for 
a large nucleus, it is not related to the large-iVc limit. We need to invoke the large-iVg 
approximation only for the contribution from the four-body S'-matrix 

(e2|S[6, |e2) = S[6, S2] = S[6, ia(r)] S[b, '^a{r')] S[6, |a(s)] S[b, |a(s + r - r')] . (95) 

Our aim is to find a momentum space representation in terms of the nuclear unintegrated 
glue. The latter is given by the pertinent function (f)g defined through 

1-S[b,^(r{r)] = J d^K^g{b,K,){l-exp[iKr]) . (96) 

Notice a subtlety: 

S[b, lair)] = exp[-la(r)T(6)] = exp[-^a,,(r)T(6)], (97) 

where we used the relationship between the dipole cross section a{r) defines for a gluon-gluon 
system and the quark-antiquark system, 

^(r) = %Ur)- (98) 

Consequently, the collective unintegrated nuclear glue (f>g{b, k) is different form the nuclear 
unintegrated glue 0(6, k) that enters deep inelastic scattering as well as diffractive quark- 
antiquark jet production and has been introduced in 0, 0]. It should not be mixed up with 
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yet another quantity, the unintegrated glue defined through the color-singlet gluon-gluon 
probe, (f)gg{K), which is defined through (for more discussion see Appendix D) 

1 - S[&, (T(r)] = J d^K(j)gg{b, Av)(l - exp[iKr]) . (99) 

If we denote by 

ao = ^ / d'Kf{K) (100) 

the dipole cross section for a large gluon-gluon dipole, then 0^^, 0^ are related as 

(f)gg{b, k) = 2S[b, ^ao](pgib, k) + (^(pg (g) 0^) (b, k) 
(03 ® 0g) (b, i^) = J d''qcj)g{b, K - q)(f)gib, q) . (101) 

Another useful quantity is 

^g{b, 1^) = S[b, |ao](5(2)(A^) + 03(b, k) , (102) 

in terms of which, 

S[b,a{r)] = J d^K<^g{b,K) exp[iKr] . (103) 

For later applications we shall also use the collective glue for a slice < /? < 1 of a 
nucleus 

S[6, Pa{r)] = J d^K $g(/3; b, k) expfi^r] . (104) 



It has the convolution property [^g{[3i) ® $g(/92) j (6, k) = $g(/5i + (32;b,K) and we note, 
that eq. (jlOU) amounts to ^gg{b, n) = $3 (2, 6, n). 

Intranuclear attenuation of dipoles becomes manifest in the distorted wave functions, 
defined in dipole and momentum space, respectively, as 

^{P;z,r) = S[b,/3a{r)]^{z,r) 

^'(/?;2,p) = Jd'^rexp[-ipr]'^{P;z,r) = J d'^K^g{b,p - k)'^{z,p) . (105) 

Now, using our master formula, the inclusive dijet cross section for gluons in the octet final 
state unfolds as 



dcr{g* {gig2}8A+8s) 1 



d^bdzd^pd^A {2ny 



J d'^sd'^rd'^r' exp[—iAs] exp[— zp(r — r')] 
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{vl>(l; tW{\- z, r') S[6, \a{s + r- r')]S[b, \a{s)] 



+^{z, r)^*{z, r') S'[b, \a{s + z{r - r'))] 



-^(1; z, r)m\z, r') S[b, \a{s - zr' + r)] S[6, \a{s - zr')] 
-m{z, r)^*(l; z, r') S[b, \a{s - r' + zr)] S[6, \a{s + ^r)]} 




X 



Ki) - ^!{z,p- ZHi - Zhi2)\^ + p - K2) - ^(^.P 



— 2Ki — 



(106) 



Here we presented the result in a manifestly ni K2 symmetric form. This result fully 
conforms with the concept of universality classes introduced in Refs. P, 0] and must be 
compared to the dijet spectra of other reactions in which the dijets are produced in the same 
color representation as the incident parton: excitation of color-octet quark-antiquark dijet 
or of open heavy flavor from gluons, g —>■ {QQ}s, and excitation of color-triplet quark-gluon 
dijets from quarks, q — > {qg}3- The incoherent distortion factor ^g{b, Ki)^g{b, K2) in the 
integrand is the same as in another gluon induced reaction g —>■ QQ. It must be compared 
to the distortion factor A) in the case of g — {qg}3- Following 0,0] we note, that (i) 
to the considered large- Ai'c approximation the above $^(6, k) equals the nuclear collective 
glue k) defined via the quark-antiquark dipoles, (ii) the incident gluon behaves as a 
pair of color-uncorrelated quark and antiquark propagating at the same impact parameter 
and (iii) $c,(&, K2)^g{b, ki) can be considered as a product of uncorrelated distortion factors 
of the quark and antiquark. As in the case of g — {qg}3, we can treat 



which contains the coUinear singularity of the g ^ gg splitting, as an intranuclear-distorted 
hard fragmentation function of the quasielastically scattered gluon. In close similarity to 
the excitation q {qg}^-, one of the wave functions which enter the fragmentation function 
is coherently distorted over the whole thickness of the nucleus. The only difference is, that 
with the incident gluon, the argument of the coherently distorted wave function is p — Hi 
versus p — A in the case of the incident quark in q {qg}^- 
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B. Coherent difFractive dijets and the rapidity gap survival 

The coherent diffractive back-to-back dijets form a universahty class of their own 
They can readily be isolated from the generic octet final states upon the application of the 
expansion (|1(J2|) : 

The incident gluon has a net color charge, and in close similarity do diffractive excitation 
of color-triplet quark-gluon dijets from quarks, q — > {qg}-^-, the coherent diffractive cross 
section is suppressed by a nuclear attenuation factor. For the color-octet incident gluon 
the attenuation is stronger than the one for the color-triplet quark, S\p^a^^\qg^ the two 
attenuation factors are related by 

S[b,a\^\,= iS[h,af\,)''^"'r (108) 



Here we notice that the nuclear attenuation of coherent diffraction can be identified with 
Bjorken's rapidity gap survival probability To this end, the relationship (jlOHj) . in 

conjunction with an absence of similar nuclear attenuation for coherent quark- ant iquark 
dijets in DIS, implies a strong breaking of diffractive factorization: the pattern of nuclear 
suppression of diffraction depends strongly on the hard subprocess and one can not treat 
diffractive production off nuclei as a hard interaction with partons of a universal nuclear 
pomeron. A more detailed discussion of this issue and possible implications for diffractive 
processes at proton- (anti)proton colliders will be reported elsewhere. 

C. Excitation of gluon-gluon dijets in symmetric octet states 

It is interesting to evaluate the contribution of symmetric octets alone. Here the relevant 
excitation operator is 

-^((e^l - (e3|)S[6,sW](|e2) + \e,))^ = i(S[6, S,] - S[6, 
= is[b, |Si](S[6, \cr{s)]S[h, \a{s + r- r')] 

-S[h,\a{s + r)]S[b,\a{s-r')]], (109) 
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which results in the cross section 



d^bdzd^pd^A 4(27r)2 



J d"^ Kid"^ K,2S {A. — Ki — H2) 



x$3(6, Ki)$3(6, K2)|*(l; z,p- Ki) - ^(1; z,p- K2)r 



4(27r) 
1 

2(2^) 



]-- J d^K(l)g{b, K)0g(6,A-At)|*(l;z,p-K)-^(l;z,p- A + k) 



+ 7T7T^ ■ S[b, |^o]..0p(6, A)|vl/(1; - 11/(1; - A)|2. (110) 



Of course, the coherent diffractive excitation of the symmetric octet state is not allowed. 



D. Excitation of dijets in higher multiplets: decuplets, 27-plet, R-j 

In order to isolate excitation of higher multiplets we have to evaluate the matrix element 
--T rf/5-^(e4|S[6,(l-/3)S(°)]cI.S[6,/?S(°)](|e2) + |e3)) 

T /'rf/5(e4|S[6,(l-/3)S(°)]|e4) {e,\Cj\e2) (e2|S[6, /?E(°)]|e2) 



2^2 Jq 

= \T{b)n{s, r, r') /' d/3S[b, (1 - P)^4S[b, (3^2] 
4 Jo 

^^T{b)nis,r,r')DAib,s,r,r') (111) 

Let us concentrate on the nuclear distortion factor 

DA{b,s,r,r') = /' S[b, (1 - /5)S4]S[6, /3S2]. (112) 
Jo 

It is of precisely the same form as excitation of color octet quark- ant iquark dipoles in 
DIS or excitation of sextet and 15-plet quark-gluon dipoles in quark-nucleus collisions. 
Here S2 describes the initial state interactions (ISl) in the slice [0, (3] of the nucleus, 
whereas S4 describes the final state interactions (FSl) in the slice 1]. Now use that 
fact that E2 = i(Si + S4) and S[6,/5E2] = S[6, i/3E4]S[b, i/3Ei]. Here we identify 
S[b, = S[b, j3^a{r)]S[b, /3|cr(r')] with the coherent distortions of the gluon-gluon dipole 

wave function in the shce [0,/5], whereas S[6, 1/5S4] will give incoherent ISl effects in the 
slice [0,(3]. 

The distortion factor pi2|l is of precisely the same form as in the excitation of color octet 



quark-antiquark dipoles in DIS j2| or in the excitation of sextet and 15-plet quark-gluon 
dipoles in quark- nucleus collisions |^]. The only difference is that both the ISl and FSl 
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cross section operators are proportional to one and the same S4, so that ISI and FSI can be 
lumped together, 

S[6, (1 - /3)S4]S[6, i/?E4] = S[b, i(2 - P)J:,] = S[b, (1 - P) + p)^,]. (113) 

Of course, in the considered large- A^c approximation we have C2[27] = 2Ca, hereafter we keep 
C2[27] on purpose as a reminder that collective nuclear glue is a density matrix in the color 
space, for which reason the /3-dependence of the nontrivial effective shce P + {1 — P)C2[27] / C a 
is controlled by the color properties of the initial and final state partons, for a related 
discussion see Ref. . The same comment is relevant to the case of decuplet dijets to be 
considered in the next section. This gives the Fourier representation 

J (fKidpK2 exp[is(K;i + 1^2)] exp[zK;2(T' — f')] 

= £ dpS[b,Pla{r)]S[b,Pl(r{r')] j £nid^K2CpK^cpK,^ 

X exp[is(K;i + ^2 + + i^a)] exp[z(K;2 + K,A){r - r')] 
X P);b,n.,Yg{P,b,K,) 

X $J^^(l-/?);b,/^4) &,/^2). (114) 

The second form emphasizes the distinction between the ISI interactions (the transverse mo- 
menta 1^1,2) and FSI (the transverse momenta K3,4) which is obscured in the first, convoluted, 
form of the distortion factor. 

That gives rise to the cross section 

d(y{g* ^ {^i^2}io+To+27+i;7) _ 1 



X 



X 



d?bdzd?pd?/^ 8(27r)2 Jq 

X j K/^d^ K^d^ nd"^ Kid^ K 5^'^\A. — k ~ ki — K2 — — K4) 

x/('^)'^'.(%P(l - b, f^3)^M b, ,,i)$,(^^(l - Py, b, b, K2) 

x||\E'(/5; z,p - ^2 - ^^4) - ^(/3; z,p- K,2- K,^- k)|^ 
+ !$(/?; z,p-K,i- K3) - ^{P; z,p- Ki- k)P} 

= I' dp J d''K,d^K,d^K2f{K,)5^^\A -^-^,-^,) 
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X % (^(1 b, (^(1 b, K,) 

X [\^{/3; z,p- Ki) - ^{f3; z,p-Ki- k) 1^ 

+ Z,p- At2) - - - A^)r} . (115) 

Again, we have a full agreement with the concept of the universality classes introduced 
in 0, 0|: the pattern of coherent distortions of the wave function and of the incoherent 
ISI and FSI distortions repeats that of other processes with excitation of dijets in color 
representations with the dimension higher by the factor oc than the dimension of the 
color representation of the incident parton. 



E. Excitation of decuplet dijets 

The case of the decuplet dijets is exceptional because they are not excited off free nucleons 
via lowest order one-gluon exchange. A new feature is that intranuclear rescattering makes 
the production of gluon dijets in the decuplets possible and it is interesting to look at their 
contribution separately. Using |1010) + |1010) = |e4) + jes), the production of decuplets is 
induced by the excitation operator 

-^^T(6)((e4| + (esl) J^' d(3S[b, (1 - .2; S[6, (le^) + les)) 

= \nbMs, r, r')D''+~'ib, s, r, r') , (116) 

with the nuclear distortion factor 

DT^{b,s,r,r') = |^'rf/?{(e4|S[b,(l-/?)E(°)]|e4)(e2|S[b,/5S(°)]|e2) 
-(e5|S[b, (1 - /3)sW]|e5)(e3|S[b,/3SW]|e3)} 
= £d(3S[b, 

x{S[6, (1 - /?)S4]S[6, i/?S4] - S[6, (1 - P)E,]S[b, i/^Ss]} . (117) 

The distinction between the ISI in the slice [0,/5] and FSI in the slice 1] is obvious. 
Repeating the analysis in the preceding Section, we readily find 

DT~\b, s, r, r') = d(3S[b, lf3a{r)]S[b, l(3a{r')] 

X {S[b, i(2 - f3)ais)]S[b, i(2 - /3)a(s + r - r')] 
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S[b, 1(2 - /3)ais + r)]S[6, i(2 - - r')]} 

-| (^exp[i/tir] — exp[z«;2^']) ( exp[— i/tir'] — exp[— iK;2r'] j | 



X 

2 

X /'rf/3S[b, |/3a(r)]S[b, \(3a{r')] 







x$,(2 - /5, b, A^i)<l>,(2 - /5, 6, K^) . (118) 

Here 

2-/5 = ^^(l-/5)+/3 (119) 

and for the sake of brevity we made an explicit use of C2[10] = 2Ca- The resulting contri- 
bution of the decuplet final states to the cross section can be cast two ways 



X 



fdl3 $,(2 - /5, 6, /^i)$3(2 - /?, b, A^2) 
{|^'(/?,2;,p-K2) -^(A^,p-«^i)P + |^(A^,P-K-K2) -^(/9,2,p-K-«^i)P 
+ |^(/?,2;,p-Av2) -^(/9,2;,p-K-K2)r + |^(/5,^,p-'«i) -^(A^,p-/«-«^i)P 

-|^(/?,2,p-K2) -^(/9,2;,P-K-Kl)|^- |^(Az,p-Ki) -^(/?,2,p-K-K2)|^} 



8(27r) 



^T(6) 1 (fn(fi^i(fK2 + Ki + K2 - A) /(At) 



X 



$,(2 - /5, 6, Ki)$3(2 - (3, b, H,2) 

Jo 



2 



(120) 



The product $^(2 - /S, b, Ki)^g{2 - /?, 6, K2) contains the term S^[6, 1(2 - /3)ao]6{K,i)6{K2). 
The similar term was the source of coherent diffraction into the antisymmetric octet. As 
one would have expected, such a coherent diffractive contribution to the decuplet cross 
section vanishes. The production of the decuplet double gluon states should lead to some 
interesting, physical consequences in the sense of the final hadron state. In Ref. j^] (see 
also 22l|2S|) it was pointed out that the process g {gg}io should lead to the production 
of baryonium states and baryon number flow over a large rapidity gap. This observation 
is based on the fact that in terms of the triplet color spinor indices the decuplet state is 
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described by the color wave function with three indices vl/*-^*^. For this reason the color 
neutralization of the decuplet state during the hadronization stage requires picking up from 
the vacuum three antiquarks (if one neglects purely gluon color neutralization which should 
be strongly suppressed as compared to the light quark mechanism due to a large effective 
gluon mass in the QCD vacuum m^-^-^ ~ R^^, where Rc ~ 0.27 fm is the gluon correlation 
radius in the vacuum In the string model js^ the hadronization proceeds through the 

breaking of the triplet color strings (three string for the decuplet state) due to Schwinger 
production of qq pairs in the color-electric field of the triplet strings [31]^. For the decuplet 
double gluon state the baryon number Nb = — 1 will be compensated by production of 
an additional baryon in the nucleus fragmentation region. In the case of the antidecuplet 
state we have Nb = 1 in the double gluon rapidity region, which can be viewed as a flow 
of valence baryon number from the nucleus region to the double gluon region. This effect 
may be important for the baryon stopping in pA and AA collisions. The corresponding 
numerical estimates will be given elsewhere. In Ref. it was pointed out that a similar 
mechanism with gluon splitting into color decuplet /antidecuplet double gluon states in the 
quark-gluon plasma produced in the initial stage of AA collisions should increase the high-p^ 
baryon production in ultrarelativistic heavy ion collisions. Numerical calculations J34| show 
that this mechanism (and the processes q — {qg}^ and q — > {qgje, discussed in p, which 
also lead to baryon production) may really play an important role in the anomalously large 
baryon/meson ratio observed experimentally at RHIC (for the recent review see j3]). 



VII. CONCLUSIONS 



We derived the nonlinear /c^-factorization for the last missing pQCD subprocess - pro- 
duction of hard gluon-gluon dijets in gluon-nucleus collisions when the nuclear coherency 
condition x < Xyi ~ 0.1 • A~^^^ holds. Although of limited importance at not so high energies 
of RHIC, this subprocess will be a dominant source of mid-rapidity and proton-hemisphere 
dijets in pA collisions at LHC The principal technical novelty is a solution of the rather 
involved seven-channel non-Abelian evolution equations for intranuclear propagation of four- 

^ In the string model for Nc = 3 baryon is usually described by the Y-configuration of three triplet strings 
connected in the so-called string junction, which plays the role of a carrier of the baryon number |32| . 
The baryonium state is a system of junction-antijunction connected by three triplet color strings. 
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gluon states. Our results for the gluon-gluon dijets in all color representations are presented 
in the form of explicit quadratures. The concept of universality classes jsf |^ is fully corrob- 
orated: The nonlinear k± factorization properties of excitation of digluons in higher color 
representations are identical to those in excitation of color-octet quark-antiquark dijets in 
DIS and quark-gluon dijets in higher color representations in qA collisions. Similar nonlinear 
k± factorization properties are exhibited by excitation of dijets in the same color represen- 
tation as the incident parton: g {gg}8, 9 {o'^js, Q ^ {o'S'js- In both g ^ gg and 
q ^ qg processes coherent diffractive excitation of incident partons with net color charge is 
suppressed by a nuclear absorption factor which can be identified with Bjorken's gap survival 
probability. The gap survival probabilities in the two cases are different. Furthermore, the 
related absorption is absent in diffractive 7* — > gg, which is indicative of a strong breaking 
of diffractive factorization. We mentioned mid-rapidity to proton hemisphere gluon-dijets 
dijets in pA collisions at LHC as a future application of the derived formalism. Still an- 
other potential application of our results for color-decuplet digluon production is a baryon 
number flow from the nucleus to large rapidity region. But, first and foremost, this work 
completes a derivation of nonlinear fc^-factorization for all pQCD processes in a nuclear 
medium and opens a way to a systematic compartive studies of high-p_|_ jet-jet and hadron- 
hadron correlations in different parts of the phase space of DIS off nuclei and hadr on- nucleus 
collisions. 

The diagonalization properties of the single-jet problem are somewhat beyond the major 
theme of this communication. Still, in view of the discussion in sec. VI. E the isolation of 
contributions from different final states is of certain interest, and we included Appendix C. 
The manifest diagonalization of the initially seven-coupled channel problem in the t-channel 
basis is, apparently, of more general interest and may find further applications in other 
problems. 

Acknowledgements: This work was partly supported by grants DFG 436 RUS 
17/101/04 and DFG RUS 17/138/05. 
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Appendix A: Useful SU{Nc) relations 



In this appendix we collect a number of identities that are helpful in the evaluation of 
the matrix elements of the dipole cross section operator. In the derivation of the projectors 
we follow closely, though in slightly different notation and normalization, Ref. 25|. Many 
useful S'f/(Ai'c)-identities can be found in Refs. [sgI. Is?!- 

If t'*,a = 1 . . . A'^ — 1 are 5'f/(Ai'c)-generators in the fundamental representation, the 
familiar /- and (i-tensors are defined through 



(A.l) 



or, 



ifabc = 2Tr([r,t^]t^) , dabc = 2Tr({r,t''}r) 



(A.2) 



The S'?7(A''c)-generators in the adjoint representation are 



{T°')bc — ifbac 



(A.3) 



so that their defining property [T"^, T^] = ifatcT'^, and the SU{Nc) transformation properties 
of (i-symbols give rise to the Jacobi identities 



'^fkam'^fmbl '^fkbrrdfmal '^fabm'^fkml 



(A.4) 



fkamdmbl dkbmfmal f abmdkml • (-^■^) 

To evaluate contractions of multiple /- and (i-symbols one makes use of the Fierz-identity 
for the fundamental generators, 

(t'^y^it-f^ =U\5'] - ^5]5t . (A.6) 



c 



They entail, that 



Tr(At")Tr(5t") = ^Tr(Afi) - ^Tr(A)Tr(5) 



Ti(At''Br) = -Tt{A)Tt{B) - -^Tt{AB) . (A.7) 



In conjunction with eq. ()A.2|) one can then obtain 
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faij fbij -^c ^ab ; 

daijdbij = — j;^ Sab , (-^-8) 

faijdbij . 

fiajfjbkfkci ^fabc j 

f f --^rl 

Jiaj Jjbk^^kci 2 <ibc 

- 4 

fiajdjbkdkci 2iV '^'^^'^ ("^■'^) 

- 12 

diajdjbkdkci 2iV ^'^^'^ ' 

fkanfnbmfmdlflck = ^ac^bd + ^ab^cd + —^{dackdkdb + dabkdkcd — dadkdkbc) 

dkandnbmdmdldlck — [OacObd + OabOcd J H [dackdkdb + dabkdkcd ) 

-^dadkdkcb- (A. 10) 

It is helpful to analyse the box and twisted-box traces of four fundamental generators, 
Rf^ = 4 Tr(t"t¥t'=) , Qed = 4 TT{n'^tY) . (A. 11) 

From eq. ()A.l|) . we obtain immediately 

J, 1 1 i 

^cd = -j;^SabScd + -{dabkdkdc + ifabkifkdc) + -{dabkfkdc + fabkdkdc) , (A. 12) 

on the other hand, 

R"^ = 4 Tr(rtVf ) = 4 Tr(t'^t''tV) , (A. 13) 

so that also 

J, 1 1 i 

Red = J^SacSbd + - {dcakdkbd + ifcakifkbd) + -{dcakfkbd + fcakdkbd) ■ (A- 14) 

We can equate the real and imaginary parts of R from eqs. ()A.12jA.14|l separately, and thus 
obtain the identities 

dabkfkdc + fabkdkdc = dcakfkbd + fcakdkbd i (A. 15) 
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2 

■J^i^ac^bd — ^ab^cd) + d^akdhbd — dabkdkdc = ifabkifkdc — ifcakifkbd 

= ifadkifkbc- (A. 16) 

For the tensor Q we have 

J, 1 1 i 

Qcd = J^^ad^bc + -^{dadkdkbc + ifadkifkbc) + Ijidadkfkbc + fadkdkbc) 

= J^i^ad^bc + '^ac'^fed " ^ab^cd) + ^^dadkdkbc + d^akdkbd — dabkdkdc) + ^^cd ) (^•1''') 

where we introduced a shorthand notation 

^^cd' = -^{^"-dkfkbc + fadkdkbcj 

= '-{fcakdkdb + d^akfkdb), (A. 18) 

and made use of eqns. ()A.15IA.16|) . 

Appendix B: Derivation of the projectors onto irreducible representations 

Our task is now to find the irreducible representations (j29p for the product of two adjoints 
and to construct the relevant projection operators. The auxihary tensors: 

^cd = 2 (^ic^fcd + SadSbc) , -4^^ = - (SacSbd - Sad^bc) (B.l) 

decompose the product representation space into its symmetric and antisymmetric parts, 
respectively: 

^ 5ac5bd = + (B.2) 

In addition, also 

[Dt\'cd = dackdkbd , [-D^jcd = dadkdkbc , [D = dabkdkcd (B-3) 

will prove helpful. All the above defined tensors S, A, D^, Dt, D^, as well as iY of eq. ()A.18|) 
are hermitian, i.e. (C^)"^ = (O^t)* = ^cd- The S'f/(Ai'c) -projectors onto the singlet as well 
as the two adjoint multiplets have up to the normalization factors manifestly the same form 
as their well-known A^^^ = 3 counterparts. 

^[1]:. = j;^^5ab5cd (B.4) 
P[8s]:d = j^^dabkdkcd = [Ds]fa (B.5) 
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project onto the symmetric singlet and octet, respectively, while 

-P[8^]cd = -J^fabkfkcd = J^ifabkifkdc (B-6) 

projects onto the antisymmetric octet. It is easily checked that they are indeed hermitian 
and satisfy the requirement 

{p[Rm = pwtpm'id = pwti (B.7) 

Using the identities ()A.9IA.10|) one finds the number of states they propagate: 

TrP[i?,] = P[R,]f, = dim[i?,] , (B.8) 

explicitly, 

TrP[l] = 1, TrP[8A] = TrP[85] = iV^ - 1 , (B.9) 

The symmetric and antisymmetric parts of our space contain, respectively 

Tr5 = \nI{NI - 1) , Tr^ = ]^{Nl - l){Nl - 2) (B.IO) 

states, so that now the problem arises to find the decomposition into irreducible represen- 
tations of the subspaces that belong to the projectors 

S^=S- P[l] - P[^s\ , A^ = A- P[8a] . (B.U) 

This is done most straightforwardly, following jsS | , by investigating the above defined tensor 
Q ()A.17jl . which takes the form 

First notice that its symmetric (antisjnnmetric) part is purely real (imaginary), 

SQS = ^eQ , AQA = i'^mQ , (B.13) 

and furthermore 

SQS = SQ = QS, AQA = AQ = QA. (B.14) 
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Now evaluate its square, - best by starting from its definition as a trace of fundamental 
generators - 

= SacSM + ^^cAfe - R^, , (B.15) 

so that 

From here we can conclude, that on the subspaces under investigation, Q"^ acts as the unit 
matrix: 

S^Q^ = Q^Si_ = Si_ , A±Q^ = Q^A± = A± , (B.17) 

therefore both subspaces decompose into orthogonal eigen-spaces belonging to eigenvalues 
±1 of the operator Q. We can then write down the projection operators 

P^^ =^it±Q)A± = ^it±iQmQ)A± (B.18) 



- (1 ± Q)S^ = -{t±^e Q)S^ . (B.19) 



To evaluate the dimensions of the associated representations, we first derive more explicit 
forms of the projectors. We start with the symmetric case. From the relations ()A.9IA.10|) 
we obtain 

2(N'^ -4) N"^ - 12 

{D^ + A)P[1] = ^^^^[1] > + A)P[85] = ^^P[Ss] , (B.20) 

and, trivially, 

{D^ + Dt)S = Du + Dt, (B.21) 

so that 

(B.22) 
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Now, 



(N^ -A)(N^ - 1) 
Tr A = , TrD„ = ^-^ ^ , (B.23) 



so that we can easily establish, that 

^ Nm.-im.+3) _ N^N.^im'S) ^^^^^ 

For iV = 3, we have 



= 27 , TrPs , 



0, (B.25) 



so that from now on we shall denote Pg^ = P[27], whereas for the other symmetric repre- 
sentation, which vanishes for iVc = 3 we shall use the notation Pg^ = PlRj]. It is interesting 
to note, that the vanishing of Ri for three colors can be related to a well-known accidental 
cancellation, namely 

P[R7]cd ^ _3 = ^{g (<^ac<56d + ^ad^bc + (^afe^cd) " [dackdkbd + dadkdkbc + dabkdkcdtj'^ , 

(B.26) 

is identically zero for SU{3) Q- 

This completes the reduction of the symmetric part, where we have 

S = P[l] + P[8s] + P[27] + P[R7] . (B.27) 

We now turn to the antisymmetric part of the product representation space, where we deal 
with two complex conjugate multiplets. Here we see, that 

QP[8a] = iQmQP[8A] = iYP[8A] = , (B.28) 

and hence 

Pf, =l{t±Q)A± = lil±iQmQ)A± 
2 



l-(A-P[8A]±iY) . (B.29) 



As Tr = 0, we have 



TtP^^ = ^ ^ , (B.30) 
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FIG. 7: The matrix element Cf of the crossing matrix. 



for SU{3) 



so that from now on 



TrP 



10, 



(B.31) 



= P[10] , PX, = P[TO] 



(B.32) 



Of course in this context it is merely a convention which multiplet we address as the decuplet 
and which as the antidecuplet. This completes our reduction of the antisymmetric part, 



A = P[8a] + P[10] + P[TO] 



(B.33) 



Appendix C: Eigenstates for the single-particle problem and the crossing matrix 

In sec. IV we computed the dipole cross section operator in the s-channel basis of 
states where gluons 1, 2 and 1', 2' respectively were in a definite 5'C/(A^c)-inultiplet. For the 
remainder of this section let us denote these states by \RR)s- 

\RR), ^ \{[9''{bi)^g\b,)]^^[g%b[)^g\b',)]_}^) 

= ^=i= P[i?]e1 \9^{bi) ® 9\b,) ® ® g'{b',)) . (C.l) 

^dim[i?] 

We mentioned, that for purposes of the single particle spectrum and total cross section, the 
dipole cross section operator is diagonalized in the i-channle basis where gluons 1, 1' and 
2,2' respectively are in definite color multiplets, we shall denote this basis by \RR)t- 

^ P[R]Z \g%b,) ® g\b,) ® ® . (C.2) 



^Jdim[R] 
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The proof, that the basis ()C.2|) indeed diagonahzes the dipole cross section operator of 
the single particle spectrum proceeds as follows. First, the dipole cross section matrix w.r.t. 
the basis ()C.2|) is obtained from the one w.r.t. ()C1|) by the swap of impact parameters 



{b,,b2,b[,b',}^{b,,b[,b2,K}. 



(C.3) 



That means, for the off-diagonal piece we obtain 



t{R'R'\a^^^\RR)t oc Q{bi, b[, 62, ^2) = (r{b[ - 62) + (r{bi - b'^) - a{bi - 62) - a{b[ - b'^) . 

(C.4) 

Now, for the single particle spectrum, one would integrate out, say P2 in the master formula 
(fT^. and in effect put 62 = b'2, but then. 



(C.5) 



Hence, for the purposes of the single particle spectrum, the off-diagonal elements of the 
dipole cross section vanish identically in the basis of states \RR)t, which is what we set out 
to prove. 

As the dipole cross section matrix in the basis ()C.2|) is obtained from the simple swap 
(1C.3j) . we can immediately give its eigenvalues Xf. 



Ai = cr(r — r') 
A2 = A3 = - 



cr(r) + cr(r') + crfr — r'] 



A4 = A5 = a(r) + 
iVc + 1 



As 
At 



iVc- Ir 



a(r) + a(r'] 



Nr. 



air I + cr r 



air — r 



(C.6) 



where, as throughout the main body of the text, r = bi — b2 , r' = b[ — b'2. The system of 
eigenvectors which belong to the A, is 

|Ai) = IA2) = |8yi8^)t, IA3) = 18585)4 

IA4) = |10T0)t , IA5) = |T010)i , lAe) = |2727)t , IA7) = \R-jR7)t (C.7) 

Clearly, once the spectrum of a matrix is known, the Sylvester formula would allow one 
to calculate any function of the matrix without knowledge of the eigenstates. In practice 
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however explicit knowledge of the latter is helpful. To obtain the color-wavefunctions of the 
states (jC.2|) . we need to establish the Fierz-type identities: 

Ptm = j2ciPsm, (C.8) 

i=l 

i.e. 

pmz = j:cipmf„ (c.9) 

i=l 

The t-channel projectors thus read, component- wise: 

P,[R]^^^Ps[R]Z. (CIO) 
The crossing matrix Cj is now obtained as (for a diagrammatic representation, see Fig.((Zj)) 

pjRMd ■ p{R^ra^ _ p{R,]z ■ pmii 
' pmf,pmii dim[i?,] • ^ ■ ' 

Apart from the complex, but hermitian structure iYs = Ps[10] — Ps[TO], explicitly 



'^O^s)cd — 2 (/cafcC^fedb + dcakfkdbj , (C.12) 

which already appeared in the decuplet projectors, the full set of color-singlet four-gluon 
states includes two more complex, but hermitian, tensor structures 

^(^i^'*)cd = -^{^fbakdkcd + dbakfkcd) i 

^iZi-^)fd=l{hakdkcd-d,akfkcd). (C.13) 

These new tensors iZ^"^^ correspond to mixed |8A8s')-states. We explicitly introduce the 
normalised states 



I(8a85)+). = y _4)(Ar2 - i) ^(^^^)c^l^"(bi) ®^7''(b2) ^/(b;) ^^'^(6^2)) , 
\{8A8sr)s = ^ ^^,_^^^^^,_^ y^(^i-^)g^l^"(bi) ® 9\b2) ® f{b[) ® g'ib',)) . 

(C.14) 

As mentioned in the main text, these states decouple in the dipole cross section operator 
(hke iY) from the states \RR)s relevant to our problem. It is straightforward to establish 
their crossing properties, namely 

zYt = zZ(-) , tZt^ = lYs , = . (C.15) 
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Using the projectors derived in appendix B we then obtain for the crossing matrix, in- 
cluding the complex tensors ()C.13|) . the result shown in eq. (C.16), and from there also 
the basis of eigenstates \RR)t displayed in eq (C.17). An SU{'i) counterpart of the crossing 
matrix (C.16) can be found, e.g. in 0|. Apparently, the crossing matrix could be used for 
an alternative derivation of the four-gluon dipole cross section matrix. 
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Appendix D: The menagerie of nuclear collective unintegrated glue 



A pertinent quantity which emerges in the description of hard processes in a nuclear 
environment is the collective nuclear unintegrated glue per unit area in the impact parameter 
space. It is not a single function which can be defined for the whole nucleus, for all hard 
processes of practical interest the description of the initial and final interactions inevitably 
calls upon a collective glue for different slices of the nucleus. Furthermore, such a collective 
glue must be regarded as a density matrix in the space of color representation, i.e., it changes 
: 'rom one reaction to another depending on color properties of the relevant pQCD subprocess 

m 

. One can trace the origin of these variations to the color-representation dependence 
of the color-dipole cross sections emerging in the description of these reactions. 

In the treatment of the nuclear structure function F2a{x, Q"^) and of the quark-antiquark 
dijets in DIS off nuclei it is advisable to use the collective glue x, k,) defined in terms of 
the amplitude of coherent diffractive quark-antiquark dijet production 0, [s^ 



Here M, U 



exp 



--a,,-(x,r)T(b) 



aqq[x, r 



) = J cfK,f{x, k)[1 — exp{iKr)] 



(D.l) 



(D.2) 



where 



and 



f{x,K) = _ J_ ,^^^^^2 



2\ _ dG{x, K^) 



d log 

is the unintegrated gluon density in the target nucleon. 

The so-defined collective nuclear glue admits a nice probabilistic expansion 



qqfi 



Here 



(D.3) 



(D.4) 



(D.5) 



w. 



(^a(&)) 



-Mb) 



1 J 



exp 



Mb) 



(D.6) 
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is a probability to find j spatially overlapping nucleons in the Lorentz-contracted ultrarela- 
tivistic nucleus, where 

Mb) = ^cr,,-oT(6), (D.7) 

is the thickness of the nucleus in terms of the number of absorption lengths for large dipoles, 
and we introduced an auxiliary infrared quantity - a dipole cross section for large quark- 
antiquark dipoles: 

croM = J d'Kf{K) . (D.8) 
The properly defined j-fold convolutions, 

— = / d Ki , 



1, (D.9) 



describe the collective unintegrated glue of j spatially overlapping nucleons in a Lorentz- 
contracted nucleus. They do not change from one reaction to another, the variations from 
k) to n) to (f)gg{b, k) BiTe fully described by the color- representation dependence 
of the overlap probabilities: 

^^(uA{b))=W,,,{^UA{b)). (D.IO) 
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